Reduced Basis representations of multi-mode black hole ringdown gravitational waves 
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We construct compact and high accuracy Reduced Basis (RB) representations of single and mul- 
tiple quasinormal modes (QNMs). The RB method determines a hierarchical and relatively small 
set of the most relevant waveforms. We find that the exponential convergence of the method allows 
for a dramatic compression of template banks used for ringdown searches. Compressing a catalog 
with a minimal match MM = 0.99, we find that the selected RB waveforms are able to represent any 
QNM, including those not in the original bank, with extremely high accuracy, typically less than 
1CP 13 . We then extend our studies to two-mode QNMs. Inclusion of a second mode is expected to 
help with detection, and might make it possible to infer details of the progenitor of the final black 
hole. We find that the number of RB waveforms needed to represent any two-mode ringdown wave- 
form with the above high accuracy is smaller than the number of metric-based, one-mode templates 
with MM = 0.99. For unconstrained two-modes, which would allow for consistency tests of General 
Relativity, our high accuracy RB has around 10 4 fewer waveforms than the number of metric-based 
templates for MM = 0.99. The number of RB elements grows only linearly with the number of 
multipole modes versus exponentially with the standard approach, resulting in very compact repre- 
sentations even for many multiple modes. The results of this paper open the possibility of searches 
of multi-mode ringdown gravitational waves. 



I. INTRODUCTION 

We are quickly approaching the era of gravitational 
wave astronomy with a world-wide network of advanced 
interferometric detectors. By 2014, all advanced LIGO 
interferometers are expected to be operational and will 
be joined by at least three other detectors, the GEO600, 
advanced Virgo, and LCGT interferometers. Detection 
rates for binary neutron star coalescences have been es- 
timated to likely be around ~ 40 yr -1 and ~ 20 yr -1 for 
binary black holes [TJ [5] . With the potential for discovery 
and rich science, it will be important to accurately detect 
these signals. The method of matched filtering, which is 
currently used for searches of known gravitational wave- 
forms in LIGO's S1-S6 runs, is the optimal method to de- 
tect signals buried in Gaussian noise [3]. While LIGO's 
noise is not Gaussian, the method, supplemented with 
additional signal-based vetoes and multi-detector coinci- 
dence requirements, is quite suitable for detection pur- 
poses, nevertheless. 

Black hole perturbation theory and numerical simu- 
lations have shown that, according to General Relativ- 
ity (GR) and many alternative theories of gravity, per- 
turbed black holes experience an exponentially decaying 
and oscillatory (or ringdown) phase in which the gravita- 
tional wave signal is dominated by a series of quasinormal 
modes (QNMs) (see [4] for a review). LIGO ringdown 
searches currently assume that the waveform is domi- 
nated by the fundamental I = m = 2 mode [5]. However, 
a single mode ringdown search for black holes of final 
mass M > 10 2 M Q can miss more than 10% of events 
in both LIGO and advanced ground-based detectors [5J. 



Furthermore, parameter estimation errors can be large 
for such single mode searches when the actual waveform 
contains a second mode 6 . Two-mode searches have 
also been proposed for consistency tests of GR and the 
no-hair theorem [BHE] and for inferring information about 
the progenitors of the final black hole formation [5]. De- 
tection of the ringdown signal from an intermediate mass 
black hole may also provide information into the forma- 
tion history of these potential gravitational wave sources 
via the black hole's spin [TO] HTj . Thus, if feasible, it is 
reasonable to expect that future ringdown searches will 
look for multi-mode signals. 

However, multi-dimensional template banks (or cata- 
logs) can result in such a large number of templates that 
multi-mode ringdown searches might not be possible in 
practice. The traditional method for building catalogs 
is to compute a metric for the parameter space that one 
uses in an algorithm to place points at a given proper 
distance away from each other and thereby guarantees a 
minimal match between a signal and template [121416] . 
Using this method for a minimal mismatch of 97%, a 
two-mode catalog suitable for testing the no-hair theo- 
rem with advanced ground-based detectors would require 
roughly 10 6 templates versus 10 3 for a similar one-mode 
search. This implies a huge increase in the computational 
cost of matched filtering searches. 

A number of reduced order modeling techniques, such 
as Proper Orthogonal Decompositions, Singular Value 
Decompositions, or Principal Component Analysis, can 
be applied to represent the original template bank by po- 
tentially much fewer elements (see, for example, [171420] ) . 

Following [21] , in this paper we show that the Reduced 
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Basis (RB) approach provides a dramatically compact 
representation of multi-mode QNM catalogs. The main 
two features responsible for such savings are (1) the expo- 
nential convergence of the training space representation 
error (the square of Eq. ^ below) with respect to the 
number of waveforms selected by the algorithm and (2) 
the linearity of multi-mode ringdown waveforms. 

We summarize some of our key results in this paper. 
Additional details and an exploration of RB represen- 
tations for ringdown searches and parameter estimation 
will be presented elsewhere. 



qi can be found in Refs. [5J |S]. For example, for the 
(2, 2, 0) and (3, 3, 0) modes considered in the following 
sections we have: 
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27TM/330 = 1-8956 - 1.3043(1 - j) 01818 
Q330 = 0.900 + 2.3430(1 - j)-°- 4810 . 



II. RINGDOWN WAVEFORMS 



III. REDUCED BASIS 



Each QNM has a characteristic complex angular fre- 
quency ui£ mn where {£, m} label the angular multi- 
pole and n labels the fundamental or overtone index. 
The real part of toemn characterizes the oscillation and 
equals 2irfi mn where fi mn is called the central frequency. 
The imaginary part is the inverse of the damping time 
Trfimn/Qimn where Qimn is referred to as the quality fac- 
tor. Thus, a single QNM waveform has the form 

htmn(t) = ^(Sll-e-'^'™/*"'' cos (2irf emn t) 
r 

.(1) 

where Ai mn {Q) is the orientation-dependent dimension- 
less amplitude, r is the distance to the source, M is 
the black hole mass, and for simplicity we have set the 
phase and arrival time to zero. We work in units where 
G = c = 1 unless otherwise stated. We will ignore the 
angular dependence in the amplitude and regard Aimn as 
constants since, as we discuss in Section [V] the linearity 
of the problem implies that the reduced bases are valid 
for all values of the amplitudes and, hence, for all angles. 

Defining hi mn {t) — for t < 0, the Fourier transform 
of (l) is 
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(2) 

For our purposes, a multi-mode signal impinging on the 
detector at time t = and its Fourier transform are then 

h(t) = ^2 h imn(t), HI) = Y2 himnif), 



with h £mn (t) and h emn (f) given by Eqs. ([1]) and (|2]), re- 
spectively. 

The central frequency and quality factor for each multi- 
pole can be related to the black hole mass M and dimen- 
sionless spin parameter j through the following fitting 
formulae: 

ftmn = — + Y^Tf — , Qimn = Ql + 92 (1 - j) q3 (3) 

where values for the (£, m, n) dependent constants fi and 



The RB approach is a framework for efficiently solv- 
ing parametrized problems, representing the solutions 
in a compact way, and predicting new ones based on 
an offline-online decomposition (see |3T] and references 
therein). In the current context, the method identi- 
fies a set of parameter points such that the associated 
waveforms constitute a nearly optimal basis for accu- 
rate "spectral" expansions of any other waveform |22) . 
If the waveforms depend smoothly on the parameters of 
the problem the expansion is expected to converge very 
quickly to the original waveform as the number of ba- 
sis elements is increased. We have found that in the case 
of gravitational waves from non-spinning compact binary 
inspirals such convergence is exponential [21 . Below we 
show that this is also the case for ringdown waveforms 
(see, for example, Fig. [I]). In fact, due to the smooth 
dependence of the theory on the source parameters we 
expect the convergence to be exponential for essentially 
all sources of gravitational waves |37j . 

In what follows we denote a generic waveform, which 
can be in the time or frequency domain, by hp and 
the relevant parameters by jl, which is in general 
multi-dimensional. In the current application, ft = 
{(fe mn ,Qimn,Ae m n)} [38]. The scheme starts with a 
training space T = {ili\f = i of P numerical values of 
the parameter fx and associated normalized waveforms 
{kpiJiLi, which we call the training space catalog. (The 
underlines denote waveforms or parameters present in the 
training space.) The natural scalar product used is di- 
rectly related to the one of Wiener filtering (see, for ex- 
ample, [53]). For two functions F,G in frequency space 
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(6) 



with S n (f) the one-sided power spectral density (PSD) 
of the detector and fhiSu frequency limits implied by 
the structure of the PSD. The scalar product of (pi) is 
related to the standard overlap integral [21] by (F,G) :— 
4Re(F, G). For dcfinitencss, in this paper we use the 
Adv LIGO PSD fit of [25]. While any PSD may be used, 
as discussed in [3T], the high accuracy of the reduced 
basis implies that it can be used with other PSD's with 
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a marginal increase in the representation error. 

The training space can be constructed by any means, 
including simple random or uniform sampling, more so- 
phisticated stochastic methods |26j . the metric approach 
[T2"] . or those of Ref. [27] j f° r example. Regardless of 
the method used to populate it, the RB formalism pro- 
duces a compact and highly accurate representation of 
the training space catalog, among other things. 

Part of the output of the algorithm is a sequential 
selection of N parameter points {/2i, fa, ■ ■ ■ , /2/v} C T 
and their associated waveforms {hp 1} hn 2 , . . . ,hp N } C 
{hp.}f =1 . The set of waveforms {hp.}f =l (or a linear 
combination of them) constitutes the reduced basis. It is 
sometimes convenient [see Eq. ([8])] , though not necessary, 
to work with an orthonormal set {ei}fL 1 instead of the 

{V.}, v ,-. 

The points {//j}^. 2 are selected by a greedy algorithm 
|28| and are sometimes referred to as greedy points. A 
choice for the first one, pi (called the seed), is needed to 
initialize the algorithm. The seed can be chosen arbitrar- 
ily, and since the greedy scheme is a global optimization 
method it is insensitive to its choice; we explicitly illus- 
trate this in Section IV B| (see Fig. [I]). The algorithm 
ends when a prescribed tolerance error, defined below in 
Eq. in the representation of the training space cat- 
alog is reached. For the salient features of the greedy 
algorithm see Alg. [TJ 



where we have implicitly introduced the projection coef- 
ficients otij . Another output of the algorithm is precisely 
the set of a^- projection coefficients for waveforms in the 
training space catalog, whereas coefficients for any other 
waveforms may be easily computed. 

The greedy algorithm sequentially selects N = N(e) 
greedy points and associated waveforms until the max- 
imum error to represent every element in the training 
space catalog reaches the specified tolerance e, 



e := 



max 

fieT ~± 



(9) 



where the norm || • || is the one induced by the scalar prod- 
uct jSJ. As discussed in [21], in the limit of sufficiently 
dense training spaces the square of e is comparable to the 
minimal match (MM) through [35] 



1 - MM as P -> oo. 



(10) 



Therefore, we refer to e 2 as the training space represen- 
tation error |40| . The exponential convergence of the RB 
method implies that e can be made arbitrarily small with 
a reduced number N of basis elements (thus the Reduced 
Basis denomination) with N < P and, in many cases, 
N < P. The quantity 



C r := P/N 



(11) 



Algorithm 1 Brief description of the Greedy Algorithm 



1: Input: {p4 ,hjx.}( =1 , e 



Seed choice (arbitrary): jli 
RB = {hp,} 
i = 1 and e\ = 1 
while 6j > e do 

i = i + l 
€i = max^r 



p (i~i)ha 



kp\ 



argmax ji j er ||P( i _i)/i^ - h p \ 



Mi 

RB = RB U hjtt 
end while 



11: Output: RB and greedy points 



The reduced basis is used to approximate other wave- 
forms, whether they were in the training space catalog 
or not, through linear combinations that represent an or- 
thogonal projection (with respect to the scalar product 
(fo)) onto its span, 



hp 



P N hp + 5hp . 



(7) 



The RB approximation is PNhp and satisfies 
(Ppfhp,Shp) = 0, by construction. If using an or- 
thonormal basis {ei}f =1 then 



N N 



1=1 



i=l 



is called the compression ratio [29]. Unless otherwise 
noted, all the results that we quote are for e 2 = 10 -12 in 
order to avoid roundoff artifacts that gradually appear 
(otherwise we could choose e 2 around double precision 
roundoff, e 2 ~ 10 -14 ). Therefore, any waveform hp in 
the training space catalog equals its RB approximation 
PNhp to this level of accuracy. 

Another error of interest is associated with the RB rep- 
resentation of any waveform hp, not necessarily present 
in the training space catalog. We call this error the wave- 
form representation error, S(ft), and define it as, 



5(jl) :=\\6hp\\ 2 = l-Re(hp,P N hp) 



(12) 



(8) § 



We quantify this error through Monte Carlo simulations, 
where we randomly sample fl and compute 5(p). Re- 
markably, as discussed in Sections |IV C| and |V A| a train- 
ing space with a modest, finite minimal match will pro- 
duce a reduced basis that represents the whole space of 
QNM waveforms with extremely high accuracy. 

The greedy points are hierarchical (i.e. nested) imply- 
ing that for extra/less accuracy in the RB representation 
points are simply added/removed. Each greedy sweep is 
embarrassingly parallel and the computational complex- 
ity in going from i to (i + 1) basis waveforms is indepen- 
dent of i. In other words, the total cost of building a RB 
with N waveforms is linear in N. These features allow 
the framework to handle large training spaces (cf. Table 
I]). Another salient aspect of the RB approach is that it 
allows one to efficiently identify on the fly those param- 
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eters most relevant for numerical simulations, without 
any a-priori knowledge of the solutions (see, for example, 
[3D]). 

In the simplest, idealized conceptual matched filtering 
search one would compute the overlap between the signal 
s and every member of the training space through its RB 
representation 



1 F 



N 



(*, PNh n . ) = 4 Re ( s > ei>«« + O (e 2 ) 



(13) 



Since the au have been precomputed offline, the filter- 
ing now involves computing significantly fewer integrals 
if TV <C P, as is usually the case. 

This idealized picture gets considerably more involved 
when including external parameters. The standard ap- 
proach for maximizing the filter output over these param- 
eters involves an extremely efficient Fast Fourier Trans- 
form (FFT). However, similar to a Singular Value De- 
composition reconstruction, the RB reconstruction com- 
putational cost of the same direct FFT strategy would 
be very high and offset the savings in the reduced order 
modeling. This is discussed in detail in Ref. 3T| so we 
do not repeat the analysis here since the issue is exactly 
the same. Strategies such as those proposed in J3TJ [32] 
should be investigated in detail to evaluate the efficiency 
of RB in actual searches. 

If one attempted a matched filter search with a RB cat- 
alog by filtering each basis function against the data and 
maximizing over arbitrary linear combinations of these 
filter outputs, one might encounter high false alarm rates. 
Instead, it is important to allow only linear combinations 
that correspond to physical waveforms; see |21j for more 
details. 



IV. ONE MODE RINGDOWN CATALOGS 
A. Metric-based catalogs 

Searches for gravitational waves from perturbed black 
holes have so far used template banks with only the fun- 
damental (£,m,n) = (2,2,0) quasinormal mode (see [5] 
and references therein). The ringdown analysis pipeline 
constructs a 2-dimensional lattice template bank in the 
(/, Q) space where the mismatch between two templates 
differing in ringdown frequency by df and in quality fac- 
tor by dQ is given by the metric [S] (T3J [33] 
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FIG. 1: Representation error as a function of the number of 
reduced basis waveforms for a single mode catalog and all 
possible choices of seeds (see Section IV B |. The dark line 



shows the average and the shaded area the maximum disper- 
sion around it. Clearly, the method is robust and does not 
require any fine tuning. Notice that the exponential rate of 
convergence of the error is present from the outset. 



and dQ 2 . Starting with the smallest central frequency 
and quality factor, the algorithm places templates along 
Q in steps of 



dQ = 



ds eS Q (1 + 4Q 2 ) 



and along 



a/3 + 16Q 4 
eflog(/) in steps of 



ds eS 



(15) 



(16) 



where ds c g = Ay/ (1 — MM) and MM is the specified min- 
imal match between template and signal. An example of 
a ringdown template bank for 



MM = 0.99, 

/220 G [10,4000] Hz. 



(17) 

? 220 G [2.1187,20] (18) 



is shown in Fig. [2] as the gray points. This catalog has 
a total of 2,213 templates. From Eqs. Q these par- 
ticular parameter choices correspond to the spin of the 
black hole j in the range [0,0.9947] and mass M in the 
range [2.9744, 3025.7] M . While the (/ 22 o,Q22o) param- 
eter space is a rectangle by construction, the correspond- 
ing shape in (M,j) space is a "warped rectangle" and so 
our stated ranges are not inclusive. Except when other- 
wise noted, the values and ranges in (17) and (18) are 
our default throughout the remainder. 



B. One- mode Reduced Basis 



which assumes white noise. Currently, the population 
algorithm uses only the diagonal terms containing df 2 



We begin by using the single-mode QNM metric tem- 
plate bank of Eqs. ( |17|18| and Fig. |] as our default 
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FIG. 2: One-mode metric placement catalog (gray dots), 
which is taken to be the training space, and the subset of 
points selected by the greedy algorithm (red bars) fo r the 
(2,2,0) mode. The settings are discussed in Sections IV A 
and lTVFll 



training space for building a reduced basis. In order to 
explicitly show the robustness of the method we use all 
2,213 possible elements of the training space as a seed 
and run the algorithm 2,213 times, once for each possi- 
ble seed. The results are summarized in Fig. [T] which 
clearly shows that the method is robust and its accuracy 
and exponential convergence rate do not depend on any 
fine tuning of the seed. This is actually expected as a 
consequence of the greedy algorithm being a global op- 
timization method. As a result, in practice one chooses 
a single, arbitrary value of the seed. Notice also from 
Fig. [I] that the exponential convergence of the error is 
not asymptotical (i.e. for large number of RB waveforms) 
but is present from the outset. 

Fig. [2] shows the points in the training space and the 
subset selected by the greedy algorithm for our default 
training space representation error of e 2 = 10~ 12 . The se- 
lected points are essentially those with the largest quality 
factor Q (corresponding to the slowest decaying modes) 
and a very few extra ones with lower Q and with large and 
small central frequencies. The figure illustrates the global 
nature of the algorithm and the redundancy present in 
sampling via a local criteria. 

As the maximum quality factor in the training space 
is increased, we find that the number of RB waveforms 
grows linearly with the corresponding number of tem- 
plates in the training space, leading to a roughly con- 
stant compression ratio [as defined by Eq. (Ill] for any 
chosen e. The qualitative behavior of selected points re- 
mains the same, with the vertical right red line in Fig. [2] 
shifting to higher Q values and the reduced basis asymp- 
totically resembling a Fourier representation with a few 
fast decaying waveforms to account for the damping of 
QNMs. 

Columns 2 and 3 in Table U show the number of RB 
waveforms needed to represent training space catalogs 
with different minimal matches. The first element pop- 
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TABLE I: Number of RB waveforms (-/Vrb) needed to repre- 
sent l-mode training spaces for the (£, m, n) = (2, 2, 0) and 
(3,3,0) QNMs with different minimal matches MM. The 
training space representation error is taken to be e 2 = 10~~ 12 . 
The number of metric-based templates scales with MM as 
Nmotric oc (1 - MM) -1 for the 2-dimensional, l-mode QNM 
catalog [12] . 



ulated by the metric placement algorithm is used as the 
seed for the greedy algorithm but the results are insen- 
sitive to any other choice (cf. Fig. [I] and related discus- 
sion). 



It might appear from Tablc|T]that for MM = 0.97-0.99 
and for a high-accuracy training space representation er- 
ror of e 2 = 10~ 12 the compression ratio should be modest 
(C r ~ 2 — 4) for a single mode catalog and would only be 
significant for larger MMs and/or smaller RB representa- 
tion errors. However, as discussed below in Section [lV C| 
it turns out that the reduced basis represents any wave- 
form in the given ranges of central frequency and quality 
factor with extremely high accuracy, leading to huge com- 
pression factors even when the reduced basis is built from 
catalogs with comparatively coarse minimal matches. 

As in the case of inspiral waveforms |21) . we find that 
for any finite range of parameters (e.g., of the central 
frequency and quality factor) we can represent the whole 
continuum of waveforms within any given training space 
representation error e 2 by a finite number of RB wave- 
forms. One way of showing this (as was done in [21]) 
is by explicitly computing the number of RB waveforms 
needed to represent training spaces built from different 
values of MM. An extrapolation to the limit that the 
training space becomes the continuum space of wave- 
forms (i.e., MM — > 1) turns out to asymptote to a finite 
number of RB waveforms (see also [31] where a similar 
result is found for MM < 0.99 in the context of singular- 
value decompositions for inspiral waveforms). Another 
way, which actually has broader implications and has not 
been demonstrated before, is discussed next. 



C. Reduced Basis representation accuracy of 
arbitrary one-mode waveforms 

Every waveform from the training space catalog is rep- 
resented with an arbitrarily small error (e 2 < 10~ 12 in all 



(> 



of our cases) by a relatively small number of RB elements. 
We now discuss the accuracy of the basis in representing 
waveforms that are not necessarily part of the training 
space. 

Starting from a catalog with a given MM to build 
the basis, the worst possible scenario would be a rep- 
resentation error for some particular waveform hp of 
1 1 <57t.^2 1 1 2 ~ (1 — MM). However, one expects the accu- 
racy to be much better since the RB framework exploits 
the global structure of the template space. As we show 
below, the accuracy is indeed orders of magnitude better 
than the pessimistic 1 — MM upper bound. 

We first built a RB representation of our default train- 
ing space catalog as specified by Eqs. (17) and (18). 
We then randomly sampled more than 6 x 10 s single- 
mode QNM waveforms drawn uniformly from the same 
ranges of central frequency and quality factor (18). We 



finally evaluated the waveform representation error ( 12 ) 



for each sample, which gives the fractional signal-to-noise 
loss from approximating a waveform by its RB represen- 
tation. The results of these Monte Carlo computations 
are shown in Fig. [3j The error in representing any of 
the randomly chosen waveforms with its RB representa- 
tion is found to be smaller than 9 x 1CP 10 . In order to 
pinpoint this upper bound we refined the Monte Carlo 
simulations to draw samples from the regions of parame- 
ter space with the largest waveform representation errors 
(i.e., the areas shaded red in Fig. [3]b) . 

The average waveform representation error for a Monte 
Carlo simulation with 10 7 points, corresponding to Fig.|3j 
is approximately 4.51 x 10~ 13 , which is less than the 
training space representation error of e 2 = 10~ 12 , and 
the most frequent value is ps 2.6 x 10 -14 . 

In other words, in all of our Monte Carlo simulations, 
any waveform, not just those present in the training space 
catalog, is found to be represented by the RB with ex- 
tremely high accuracy. 

A simple extrapolation from Table|l]shows that ~ 10 10 
metric templates would be needed to achieve a MM com- 
parable to the above waveform representation error of 
r~> 1CP 9 for our reduced basis with 505 elements, imply- 
ing an effective compression ratio of C r ~ 10 7 . Further- 
more, the fact that our simulations strongly indicate that 
the maximum waveform representation error is strictly 
bounded (by 9 x 10~ 10 in the settings discussed) would 
imply a formally infinite compression ratio. This again 
reflects another facet of being able to represent the whole 
spectrum of waveforms with a finite number of basis ele- 
ments. The key distinction here, compared to the discus- 



sion at the end of the previous subsection ( IV B ) , is that 



one does not need to explicitly compute a set of bases 
with increasingly larger number of training space points 
to reach a target accuracy limit. This would allow for 
the a priori construction of high accuracy reduced bases 
for sources with larger number of parameters, such as 
precessing binary inspirals, at a fraction of the projected 
computational cost. 

We have done exhaustive additional tests to support 



the conclusion that the representation error is extremely 
small for all waveforms with parameter values in the set 
ranges [e.g., as in Eq. (18)]. For example, since a ran- 
dom, uniform sampling of (/220,Q22o) leads to a non- 
uniform sampling of (M,j) [see Eq. Q] we have per- 
formed Monte Carlo simulations with random, uniform 
sampling of mass and spin such that the correspond- 
ing central frequencies and quality factors are within the 
ranges of Eq. ( 18 ). In addition, instead of selecting wave- 
forms at random, we have used all those from a metric- 
based catalog with a minimal match of about 0.9999, cor- 
responding to more than 10 5 templates. We represented 
the waveforms using the RB generated by our default 
training space [from Eqs. (17) and (18)] and computed 
the corresponding waveform representation errors. In all 
cases the errors were found to be bounded by the previ- 
ously quoted value of 9 x 10~ 10 . 

As discussed in the next section, the high accuracy 
associated with representing arbitrary waveforms by a 
RB carries over to the case of multiple QNMs. 



V. MULTI-MODE RINGDOWN CATALOGS 

In this section we discuss how the RB approach can 
be used to efficiently build, compress, and represent the 
space of multiple QNM waveforms. We present two kinds 
of multi-QNM reduced basis. The first is for a two-mode 
ringdown waveform wherein the two central frequencies 
and quality factors are related to each other through a 
constraining relation provided by GR for distorted black 
holes. The second is for a ringdown waveform consist- 
ing of p modes where the central frequencies and quality 
factors are unconstrained by any model or theory. In 
particular, this latter case is amenable for providing a 
test of the no-hair theorem [6] , through a test of the pre- 
cise relationship between the multiple ringdown modes 
as predicted by GR. 



A. Constrained two-mode Reduced Basis 

Here we consider two-mode ringdown waveforms with 
(I, m, n) — (2, 2, 0) and (3, 3, 0), which are of the form 



h = C [(1 - A)h 220 + Ah 330 ] 



(19) 



where A G [0,1] and C is fixed by the normalization con- 
dition (h, h) = 1. 

The parameters for this two-mode QNM are 
(/220, Q220), (/330,Q33o), and the relative amplitude pa- 
rameter A. However, GR provides a relation between 
{flmm Qimn) and the mass and spin of a perturbed black 
hole. From Eqs. ([3| or, specifically for the two-mode case 
here, Eqs. Q and ^ it follows that (/33o,Q33o) are re- 
lated to (/2207 Q220) so that there are only three indepen- 
dent parameters. In this sense, GR constrains the oth- 
erwise 5-dimensional parameter space to a 3-dimensional 
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FIG. 3: Waveform representation error of randomly selected 
QNMs that are not necessarily present in the training space, 
(a) Distributions of waveform representation errors for 1- 
mode (gray) and 2-mode (white) using reduced bases built 
from a training space with MM = 0.99. In all cases, the error 
is found to be < 9 x 1CT 10 for sin gle mode s and < 9.8 x 10~ 10 
for two modes. See Section IV C and V A for details, (b) The 
waveform representation error [defined in ( 12 1] as a function 
of the randomly chosen values of (/220, Q220) for the 1-mode 
case with the same setup as in (a). The color map indicates 
the value of the error on a log 10 scale. The worst errors oc- 
cur for the largest Q220 values and are mostly independent 
of /220- The crosses indicate the parameters selected by the 
greedy algorithm for this reduced basis. 



one, namely, {/220, Q220, A} (or, equivalently, {j, M, A}). 

For a given MM and for A G [0, 1] we build a train- 
ing space for the (2,2,0) mod e usin g the metric ap- 

Next, using 



B2- 



proach as described in Section IV A 
the fitting formulae given by Eq. (|4j), the values of mass 
and spin for each template are determined. Then, using 
Eq. ([5]), the corresponding values of (/330, Q330) a re com- 
puted and paired to the starting (/220, <?22o)- This proce- 
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FIG. 4: Values for the relative amplitude parameter A in the 
two-mode constrained waveforms [Eq. ( 19 1] selected by the 



greedy algorithm for each reduced basis waveform (indexed on 
the horizontal axis). Here there are ua = 1,000 samples for 
A e [0, 1] but essentially either the (2, 2, 0) mode (.4 = 0) or 
the (3, 3, 0) mode (A = 1) are selected, with a higher density 
of points associated with the latter due to its larger maximum 
quality factor for any given black hole spin parameter value. 
About 3% of the selected amplitudes are equal to zero (not 
just nearly so). All of the remaining points selected are either 
equal to 1 (about 79%) or between A = 0.93 and 1 (about 
18%). 



dure is repeated for each template in the (2, 2, 0) catalog 
and subsequently populates the (/33o,<333o) plane with 
templates that are "inherited" from the (2, 2, 0) training 
space. The relative amplitude parameter A is sampled 
with riA equally spaced points. In this section and in 
Table [I] for the (3,3,0) mode column, MM refers to the 
minimal match of the (2, 2, 0) starting catalog. 



1. Results 

In all of our numerical simulations we have found that, 
for any fixed training space representation error e 2 , the 
number of reduced basis waveforms is independent of 714 
for 714 > 2. This is not too surprising, given the linearity 
of the two- mode waveform (19); see also the discussion 
in Section |VB| A related observation is that, as shown 
in Fig. |4j the greedy algorithm essentially selects either 
the (2, 2, 0) mode (A = 0) or the (3, 3, 0) one (A = 1) as 
opposed to combinations of them. As a result, the dis- 
tribution of parameter points selected by the algorithm 
closely resembles the union of a subset of points for two 
individual, one-mode reduced bases. This is illustrated 
in Figure [5j which should be compared to Figure [2] for a 
one-mode (2, 2, 0) reduced basis. 

The number of constrained, two-mode RB waveforms 
tVrb for different values of MM is shown in Table |H| in the 
column labeled "2-mode, GR." Notice that these values 
are only marginally larger than those corresponding to a 
one-mode (3, 3, 0) basis in Table |lj To understand why, 
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FIG. 5: The metric-based training space (points) for the same 
constrained, two-mode catalog of Fig. [4] and the subset of 
parameter values selected by the greedy algorithm (bars) for 



a training space representation error of e =10 



Red: 



(2, 2, 0) mode. Blue: (3, 3, 0) mode. The distribution of points 
closely resembles a subset of the two one-mode greedy points 
due to the linearity of the problem and to the "decoupling" 
of selected amplitude points (see Fig. p}. 



notice from Fig. [4] that the (3,3,0) mode, corresponding 
to A = 1 in ( fl9| ), is chosen much more frequently than the 
(2, 2, 0) mode. In turn, the greedy algorithm's preference 
for A = 1 can be understood by the fact that: i) for any 
maximum value of Q in the starting (2, 2, 0) catalog the 
corresponding value for the (3, 3, 0) mode for the same 
black hole spin parameter j is larger (see Eqs. Q and 
((5|)), and ii) as we have already discussed in Section V A 



the greedy algorithm for a single-mode RB mostly selects 
points with the largest value of Q in the catalog (see 
Figure [2 ) . 

Table II] also shows the approximate number of tem- 
plates iVmetric for constrained two-modes using the metric 
placement method in the column labeled "2-mode, GR." 
The compression ratio with respect to the metric number 
of templates for MM = 0.99 is C r s» 24 and dramatically 
increases as MM increases. The number of RB waveforms 
asymptotically approaches about 1,000, which suggests 
that a two-mode search with the RB approach may be 
easily feasible. 

Next, we show that the accuracy of 2-mode reduced 
bases and the associated compression ratios are much 
better, like in the 1-mode case of Section |TV C[ than ini- 
tially indicated in Table [III which are already very good. 



2. Reduced Basis representation accuracy of arbitrary 
two-mode waveforms 

In analogy with the 1-mode case, we have found that 
two-mode, constrained reduced bases built from rela- 
tively coarse training spaces turn out to represent any 



1 - MM 


2-modc, GR 


2-mode 




■*v metric 


N w 


-^metric 


Nkb 


0.03 


3.5 x 10 3 


737 


3.4 x 10 6 


1,198 


10~ 2 


1.8 x 10 4 


751 


5.3 x 10 7 


1,237 


10~ 3 


5.8 x 10 5 


958 


1.9 x 10 10 


1,495 


10- 4 


1.8 x 10 7 


1,007 


5.3 x 10 12 


1,567 


10~ 5 


5.8 x 10 8 


1,018 


1.9 x 10 15 


1,590 



TABLE II: Number of reduced basis waveforms (./Vrb) needed 
to represent 2-mode training spaces with (£, m, n) — (2, 2, 0) 
and (3, 3, 0) for different minimal matches MM. The training 
space representation error is e 2 = 10 -12 . The "2-mode, GR" 
and "2-mode" cases are discussed in Sections |V A| and |VB| 
respectively. The number of metric-based templates in those 
two cases scales with MM as AT me t r ic oc (1 — MM)~ d/2 , with 
d — 3 and d — 5, respectively [TJ]. For the "2-mode, GR" 
entries, MM refers to the minimal match associated with the 
training space for the (2, 2, 0) mode, which is used to generate 
the points in the (/33o,Q33o) plane as described in the text. 
For the "2-mode" case the training space representation error 
of 1CP 12 corresponds to each mode separately and the cor- 
responding total error is bounded from above by 4 x 10 -12 , 
as discussed in Section |VB| We are grateful to V. Cardoso 
for his help in generating the values of A r motr ic in the second 
column for the constrained, 2-mode case. 



waveform with extremely high accuracy. Since the de- 
tails follow those of Section |IV C| our description here is 
more succinct. 

In order to emphasize our finding that for constrained, 
two-mode ringdown waveforms the number of basis el- 
ements saturates at = 2 amplitude samples, we de- 
scribe waveform representation error results for a reduced 
basis built with exactly 714 = 2. That is, from here on 
the training space and, as a consequence, the resulting 
reduced basis only includes A = and A = 1 relative 
amplitudes in the waveforms. 

We have uniformly sampled the constrained, 2-modc 
parameter space (/220, Q220, A) with /220 and Q220 in 



the same range used to build the reduced basis, Eq. (18 1, 
and with A € [0,1]- The corresponding values of 
/330 and Q330 were computed as described earlier using 
the constraint from General Relativity that they corre- 
spond to the same black hole mass and spin as those 
of (/220,Q22o)- The resulting waveforms were then pro- 
jected onto the reduced basis built with a training space 
of MM = 0.99 as described above and the waveform rep- 
resentation error ( 12 ) for each sample was computed. We 



have executed many Monte Carlo simulations, with the 
largest one having 10 9 random triples (/22O1 Q220, -4)- In 
all cases the maximum waveform representation error was 
found to be below 9.8 x 10~ 10 , which is remarkably close 
to our strict numerical upper bound found for the one 
mode case described in Section IV C namely, 9 x 10 -10 . 
As an example, Fig. [3^ shows a histogram of the wave- 
form representation error for a Monte Carlo simulation 
with 10 7 points. In this case, the average error is about 
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Representation error for random two-mode waveforms 



Q220 



1 




(fimmQimn) and (p — 1) for the relative amplitudes of 
the normalized waveform. Already for p — 2 this results 
in a very large bank of roughly 10 6 — 10 7 templates when 
using the metric approach for MM = 0.97 — 0.99 (see 
Table [TT]). 

One could build a "full" reduced basis from a (3p — 1)- 
dimensional training space. However, the linearity of the 
problem and the results of Section V A suggest a much 
simpler approach. Since ringdown waveforms are a linear 
superposition of single modes [in the following the index 
7 generically labels the (£,m,n) triple], 



(20) 



i=i 



then a "simple" RB representation for ([20]) could be given 
by building a reduced basis for each 7-th mode separately 
from the 7-th training space 7} = and to define 

the representation of the multi-mode wave to be the sum 
of the individual projections 



P. 



N 



hp :— ^ AiPNjhi.rt ■ 
1=1 



(21) 



FIG. 6: Waveform representation errors, as defined by 



Eq.(12l, of constrained, two-mode waveforms for a random 
sample of (/220, Q220, A) triples. The training space used to 
build the reduced basis is relatively coarse yet the waveform 
representation error is fo und to be, in all cases, smaller than 
See Section I V A| for more details. 



9.8 x 10' 



Here Pn, denotes the standard orthogonal projection 
onto the reduced basis for the single 7-th mode built 
from a training space catalog associated with 7j. De- 
noting by e 2 the training space representation error for 
the 7-th mode, 



e 2 := max 



P: 



N 



(22) 



7.35 x 10 13 , which is less than the training space rep- 
resentation error of e 2 = 10~ 12 , and the most frequent 
value is « 2.4 x 10~ 13 . 

Fig. [6] shows the distribution of the waveform repre- 
sentation error as a function of the sample values for 
(/22O) Q220,A). As with the 1-mode case, the largest er- 
rors occur for larger Q220 values and are mostly indepen- 
dent of /220 and the relative amplitude. 



B. Unconstrained multi-mode Reduced Basis 

As mentioned in Section |TJ searches for multi-mode 
ringdown gravitational waves have been proposed as a 
consistency test of GR and, more specifically, the no- 
hair theorem. The basic idea is not to enforce the con- 
straint that the ringdown frequencies of the different 
modes should correspond to the same black hole mass 
and spin, as in the previous subsection (VA), but to 



check it a posteriori, by assuming the standard relation- 
ship in Eq. ^ between QNM frequencies and black hole 
mass and spin for each mode. 

The dimension of the parameter space for p uncon- 
strained fundamental modes [?2] is (3p — 1): 2p for 



then the total training space representation error for ( 20 ) 
is bounded by 



simple 



:= max \\h n — Pfjhj, 
peT £ —t 



1=1 



p 

£ 

j,j=i 



where the full training space T is the product of the 
individual training spaces ^j—iTi- If all the errors are 
chosen to be comparable (i.e, ej ~ e) then the bound is 
simply 



2 <- 2 2 

^simple r^i P ^ ■ 



(23) 



In most cases of practical interest for gravitational wave 
searches, including advanced and third generation earth- 
based detectors such as the Einstein telescope, p < 4 
9J. If e 2 = 10~ 12 , as in all of our cases, then the total 
training space representation error would still be very 
small, of order e 2 implo < O(l0~ n ). In fact, one could 
include up to p — 10 5 modes and still have a maximum 
total error e 2 implc < 0.01, which is comparable in terms of 
waveform representation accuracy to a minimal match of 
about 0.99. This feature is another advantage of having 
RB representations with very high accuracy, even if the 
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latter exceeds the accuracy of the detector itself or of 
the physical modeling. Such representations leave room 
for simplifications while preserving high accuracy in the 
overall approach. 

The total number N = YjJ=i of RB waveforms 
scales linearly with p. If iVj ~ J\f basis elements are 
needed to represent each mode then N ~ pAf. In con- 
trast, the number of metric-based templates is propor- 
tional to (1 — MM)^ 3 ^ 1 '/ 2 12J and increases dramati- 
cally with p for any fixed MM (see, for example, Tables 
0and|Ilj| 



Furthermore, the representation (21 1 is valid for all 
values and ranges of amplitudes Ai- In particular, N is 
independent of the range or values of the {.4/}. This 
last observation, though trivial, helps to understand the 
results of Section IV AI for constrained 2-mode waveforms: 
essentially one mode or the other is selected at each step 
by the greedy algorithm (see Fig. |4| , the number of RB 
waveforms saturates at nj^ = 2 samples (.4 = and 
A = 1), and the selected (/, Q) points (Fig. |5J) resemble 
a subset of those picked by the two 1-mode reduced bases. 

Already for the p = 2 case and template banks with 
minimal matches as low as MM = 0.97 and 0.99, this sim- 
ple RB approach for unconstrained multi-mode QNMs 
implies a compression ratio of around three and four 
orders of magnitude relative to the respective training 
spaces. In addition, as discussed in Sections |IV C| and 
V A| the resulting bases actually represent the whole 



space of waveforms with an accuracy comparable to the 
training space representation error, resulting in effec- 
tively much larger, if not formally infinite, compression 
ratios. The number of unconstrained 2-mode [(2, 2, 0) 
and (3, 3, 0)] metric-based templates for different values 
of MM and the number of RB waveforms needed to rep- 
resent them with this simple approach are given in Table 

im 



VI. FINAL REMARKS 

We have shown that the RB approach provides very 
compact and high-accuracy representations of multi- 
mode ringdown gravitational waves. For example, the 
number of RB waveforms needed to represent any two- 
mode General Relativity quasinormal mode with a rep- 
resentation error of 10~ 12 is smaller than the number 
of metric-based templates for a one-mode catalog with 
a minimal match of MM ~ 0.99. The comparison with 
the number of two-mode, metric-based templates for a 
catalog with MM — 0.97 - 0.99 designed to test the con- 
sistency of GR and the no-hair theorem is even more 
striking: the number of RB waveforms needed to repre- 
sent the continuum is 3-4 orders of magnitude smaller 
than the number of metric-placed templates for minimal 
matches of 0.97-0.99. 

We have demonstrated the robustness of the RB 
method by observing essentially identical training space 
representation errors regardless of the seed value. In the 



case of one-mode ringdown waveforms, we showed how 
the greedy algorithm preferentially selects waveforms 
with the largest quality factors to construct the basis. 
The number of RB waveforms for the (£, m, n) — (2, 2, 0) 
one-mode case needed to represent the continuum asymp- 
totes to about 600 for advanced ground-based gravita- 
tional wave detectors. We observed that the compres- 
sion ratio remains constant as a function of the maxi- 
mum quality factor. Through detailed numerical stud- 
ies we demonstrated that any waveform that falls within 
the one- and two-mode training space's ranges in cen- 
tral frequency and quality factor (and, in the case of 2- 
modes, for an arbitrary relative amplitude) can be recov- 
ered with an error no larger than s» 9 x 10~ 10 . Thus, 
we can construct high accuracy reduced basis represen- 
tations of these ringdown waveforms very easily off-line. 
We expect this result to extend to RB representations 
of spinning binary inspirals and other multi-dimensional 
systems that depend smoothly on their parameters. 

In the case of constrained two-mode ringdown wave- 
forms, we find that the RB method, in additional to its 
exponential convergence, benefits from the linearity of 
the QNMs. The number of RB waveforms is independent 
of the number of relative amplitude parameter samples 
n.4 for n_4 > 2 for any fixed training space representa- 
tion error. The greedy algorithm essentially only selects, 
for example, the (2,2,0) mode or the (3,3,0) mode. The 
number of RB waveforms for the constrained two-mode 
case to represent the continuum asymptotes to about 
1,000 for advanced ground-based detectors. Finally, in 
the case of unconstrained two-mode ringdown waveforms, 
we again exploit the linearity of the problem to propose 
the construction of a "simple" RB for each mode sepa- 
rately. Then the representation of the multi-mode wave 
is simply the sum of the individual projections. In ad- 
dition, the number of RB waveforms scales linearly with 
the number of modes indicating a very large reduction in 
the computational cost of multi-mode ringdown searches. 

The results of this paper open up the possibility of 
searches of multi-mode ringdown gravitational waves 
thereby allowing one to test the no-hair theorem, to sig- 
nificantly reduce the event loss rate of these kinds of 
signals, to improve parameter estimation, and to pos- 
sibly infer progenitors from the relative amplitudes of 
the different QNMs. These results are made freely avail- 
able ESI. 
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